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^ ■ Abstract 



We study the rare B decay B — > K*£ + £ using soft-collinear effective theory (SCET). At leading 



power in l/m^, a factorization formula is obtained valid to all orders in a s . For phenomenological 
application, we calculate the decay amplitude including order a s corrections, and resum the loga- 



rithms by evolving the matching coefficients from the hard scale 0(rrib) down to the scale \frn~b^h- 



Q_i' The branching ratio for B — > K*(. + l is uncertain due to the imprecise knowledge of the soft form 

factors C±{q 2 ) an d C||(<? 2 )- Constraining the soft form factor C-l(q 2 = 0) from data on B — ► K*j 



yields Clw = 0) = 0.32 ± 0.02. Using this input, together with the light-cone sum rules to deter- 



mine the g 2 -dependence of C±(q 2 ) an d the other soft form factor C||(<? 2 )j we estimate the partially 
integrated branching ratio in the range 1 GeV 2 < q 2 < 7 GeV 2 to be (2.921°;^) x 10~ 7 . We discuss 
how to reduce the form factor related uncertainty by combining data on B — ► p{— > itit)£v£ and 
B — > K*(—> Kir)^^ . The forward-backward asymmetry is less sensitive to the input parameters. 
In particular, for the zero-point of the forward-backward asymmetry in the standard model, we 
get q 2 , = (4-07^q'{3) GeV 2 . The scale dependence of q 2 is discussed in detail. 
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I. INTRODUCTION 



The electroweak penguin decay B — > K*£ + £~ is loop- suppressed in the Standard 
Model(SM). It may therefore provide a rigorous test of the SM and also put strong con- 
straints on the flavor physics beyond the SM. 

Though the inclusive decay B — > X s £ + £~ is better understood theoretically using the 
Operator Product Expansion, and the first direct experimental measurements of the dilepton 
invariant mass spectrum and m^-distribution are already at hand being an inclusive 

process, it is extremely difficult to be measured in a hadron machine, such as the LHC, which 
is the only collider, except for a Super-B factory, that could provide enough luminosity for the 
precise study of the decay distributions of such a rare process. In contrast, for the exclusive 
decay B — ► K*£ + £~, the difficulty lies in the imprecise knowledge of the underlying hadron 
dynamics. Experimentally, BaBar |3| and Belle 0| Collaborations have observed this rare 
decay with the branching ratios: 

f (7.8 + {i ± 1.2) x 1(T 7 (BaBar) , 

B(B -> K*£ + £~) = { y - 1 - 7 > y 1 ' (1) 

[ (16.51^ ±0.9 ±0.4) x 10~ 7 (Belle) . 

We note that the Belle measurements are approximately a factor 2 higher than the corre- 
sponding BaBar measurements. In addition, Belle has published the measurements 0, 0] of 
the so-called forward-backward asymmetry (FBA) In particular, the best-fit results by 
Belle for the Wilson coefficient ratios for negative value of Aj, 

^ = -15.332 ±1.1, 

4^ = 10.3iti ± 1.8 , (2) 

are consistent with the SM values Ag/Aj ~ —13.7 and Aiq/Aj ~ ±14.9, evaluated in the 
NLO approximation (see Table I). With more data accumulated at the current B factories, 
and especially the huge data that will be produced at the LHC, it is foreseeable that the 
dilepton invariant mass spectrum and the FBA in this channel will be measured precisely 
in several years from now, allowing a few % measurements of the Wilson coefficient ratios 
and the sign of A 7 . 



Theoretically, the exc 



see for example 



usive decay B — > K*£ + £ has been studied in a number of papers, 



From the viewpoint of hadron dynamics, the applica- 



tion of the QCD factorization approach to this channel deserves special mention, 



as we shall be comparing our phenomenological analysis with the results obtained in this 
)aper. The emergence of an effective theory, called soft-collinear effective theory (SCET) 



provides a systematic and rigorous way to deal with the perturbative 
strong interaction effects in B decays in the heavy-quark expansion. A lot of theoretical 
work has been done in SCET related to the so-called heavy-to-li ght transitions in B decays, 



in particular, a demonstration of the soft-collinear factorization 



19 



20 



22, 



21 



22, 



231 ] . a complete 



24j , and the extension 



catalogue of the various 2-body and 3-body current operators 
of SCET to two effective theories SCET/ and SCET//, with the two-step matching QCD — > 
SCET/ — > SCET// 25]. Among various phenomenological applications reported in the liter- 
ature, SCET has been used to prove the factorization of radiative B — > decays at leading 
power in l/m& and to all orders in a s Likewise, SCET, in combination with the 

heavy-hadron chiral perturbation theory, has also been used to study the forward-backward 
asymmetry in the non-resonant decay B — > Ktc£ + £~ in certain kinematic region 28]. In this 
paper, our aim is to use SCET in the decay B — > K*£ + £~ . Due to the similarity between 
B — > K* , y and B — > K*£ + £~ decays, our approach is quite similar to the earlier SCET-based 
studies j3,l3]> m particular to Ref. j^J. Moreover, an analysis of the exclusive radiative 
and semileptonic decays B — > K*j and B — > K*£ + £~ in SCET can be combined with data 
to reduce the uncertainties in the input parameters. In particular, as we show here, the lo- 
cation of the forward-backward asymmetry in B — > K*£ + £~ can be predicted more precisely 
than is the case in the existing literature. 

It is well known that, when q 2 , the momentum squared of the lepton pair, is comparable 
to Mj/i/,, the resonant charmonium contributions become very important, for which there 
is no model- independent treatment yet. Likewise, for higher g 2 -values, higher ^-resonances 
(■0', ip", ...) have to be included. Thus, in the following we will restrict ourselves to the 
region 1 GeV 2 < q 2 < 7 GeV 2 , which is dominated by the short-distance contribution. Note 
that the lower cut-off 1 GeV 2 is taken here because, as we shall see later, when q 2 is very 
small, say q 2 ~ 0(Aq CD ), the factorization of the annihilation topology breaks down. In 
this kinematic region, a factorization formula for the decay amplitude of B — > K*£ + £~, 
which holds to 0(a s ) at the leading power in l/m&, has been derived in Ref. Q] using 
the QCD factorization approach We shall derive the factorization of the decay amplitude of 
B — > K*£ + £~ in SCET, which formally coincides with the formula obtained by Beneke et 
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al. 1141. but is valid to all orders of a„: 



1. Q, 

(K:e + £-\H eff \B) = T a V)ag 2 ) + / — 0±(") / du ^(u)T a ^, u, q 2 ) , (3) 

^ Jo ^ Jo 



± 

where a =||, _L denotes the polarization of the K* meson. The functions T 1 and T 11 are per- 
turbatively calculable. ( a (q 2 ) are the soft form factors defined in SCET while 0^ and 4> a Kt , 
are the light-cone distribution amplitudes (LCDAs) for the B and K* mesons, respectively. 
Compared to the earlier results of Ref. [l4j, obtained in the QCD factorization approach, 
the main phenomenological improvement is that for the hard scattering function T n , the 
perturbative logarithms are summed from the hard scale fib ~ C?(m&) down to the interme- 



diate scale fit = vTVVi, where A h represents a typical hadronic scale. Note also that the 
definitions of the soft form factors ( a {q 2 fo r our SCET currents, defined subsequently in 
section 2, are different from those of Ref. 141 . a point to which we will return later in section 
3. Hence, the explicit expressions for T 1 , derived here and in Ref. |l4| are also different. 

This paper is organized as follows: In section II, we briefly review the basic ideas and 
notations of SCET. We then list the relevant effective operators in SCET and do the explicit 
matching calculations from QCD to SCET/ (Sec. II A) and from SCET/ to SCET// (Sec. II 
B). The matrix elements of the effective SCET operators are given in Sec. II C. At the end of 
this section, the logarithmic resummation in SCET/ is discussed. In section III, we consider 
some phenomenological aspects of the B —> K*£ + £~ decay. We first specify the input 
parameters, especially the soft form factors (± t \\(q 2 ) (Sec. Ill A), which are the cause of the 
largest theoretical uncertainty. We use the g 2 -dependence of the related QCD form factors 
in the LC-QCD sum rule approach, but fix the normalization of these soft form factors 
using constraints from data on the exclusive decays B — > K*^. In Sec. Ill B, we work 
out numerically the evolution of the B-type SCET/ matching coefficients, defined earlier 
in Sec. II. We then give the dilepton invariant mass spectrum and the forward-backward 
asymmetry in the decay B — > K*£ + £~, and compare the integrated branching ratios with 
the measurements from BaBar and Belle (Sec. Ill C). We end with a summary of our results 
in section IV and suggestions for future measurements to reduce the model dependence due 
to the form factors and other input parameters. 
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II. SCET ANALYSIS OF B — > K*£+£- 

For the b — > s transitions, the weak effective Hamiltonian can be written as 

10 



h eff ~ 



^wE^w . ( 4 ) 
* i=i 



where we have neglected the contribution proportional to V* s V u b in the penguin (loop) 
amplitudes, which is doubly Cabibbo-suppressed, and have used the unitarity of the CKM 
matrix to factorize the overall CKM-matrix element dependence. We use the operator basis 



introduced in 



mm 



(5) 



Qi = (sT A c)v_ A (cT A b) v _ A , Q 2 = (sc) v _ A (cb)v-A , 

Q 3 = 2 (sb) v _ A ^(g 7 M <?) , Q 4 = 2 {sT A b) v _ A J^i^^) , 

q q 

q 

Qe = 2 s lfl7ulp (l - 7 5 )T A 6 5>W7 P T A g) , 

q 

Qr = - 9 -^sa^(l + 75)^ , Qs = + ^ ' 

Q 9 = ^(sb) v „ A {hn) , Q w = ^(sI)v-a(W) , 

where T" 4 is the SU(3) color matrix, a em = g 2 m /47r is the fine-structure constant, and m^/i) 



is the current mass of the b quark in the MS scheme at the scale fi. 

Restricting ourselves to the kinematic region 1 GeV 2 < q 2 < 7 GeV 2 , the light K* 
meson moves fast with a large momentum of the order of mg/2, which thus can be viewed 
approximately as a collinear particle. For convenience, let us assume that the K* meson 
is moving in the direction of the light-like reference vector n, then its momentum can be 
decomposed as p^ = n-pn fl /2+p p L + n-pn fJ- /2, where n M is another light-like reference vector 
satisfying n ■ n — 2. In this light-cone frame, the collinear momentum of K* is expressed as 

p= (n-p,n-p, p±) ~ (A 2 , 1, \)m b , (6) 

with A ~ A/ nib <C 1. In addition to this collinear mode, the soft and hard-collinear modes, 
with momenta scaling as (A, A, \)rrib and (X,l, VX)m b , respectively, are also necessary to 
correctly reproduce the infrared behavior of full QCD. 



5 



SCET introduces fields for every momentum mode and we will encounter the following 
quark and gluon fields 



£ C ~A, A£~(A 2 ,1,A), e^e^-A 1 / 2 , A£ c ~ (A, 1, A 1 



/ 2 1 



A 3 / 2 , ^~(A,A,A), h ~ A 3/2 . (7) 



In the above, the symbol stands for the gluon field, h represents a heavy-quark field, 
the symbols £ and q stand for the light quark fields, and the subscripts c, s, he stand for 
collinear, soft and hard-collinear modes, respectively. Note that the momentum q of the 
lepton pair is taken as a hard collinear momentum, since in this paper we only consider 
the range 1 GeV 2 < q 2 < 7 GeV 2 . That is why an extra hard-collinear field in the 
n direction is required later. As explained in detail in Ref. [3], to construct the gauge 
invariant operators in SCET, it is more convenient to introduce the building blocks, given 
below, which are obtained by multiplying the fields by the Wilson lines which run along the 
light-ray to infinity: 

. (8) 

For example, the field Xhc is defined as 

X hc (x) = Wl c (x)£ hc (x); with W hc (x) = Pexpiig / dsn ■ A hc (x + sn) ) , (9) 



where Wh c (x) is the hard collinear Wilson line. The notations Q s and CK S are used when the 
associated soft Wilson lines are in the n-direction. For the definitions of the other fields and 
more technical details about SCET, we refer the reader to Ref. 2fJ and references therein. 

Since SCET contains two kinds of collinear fields, i.e. hard-collinear and collinear fields, 
normally an intermediate effective theory, called SCET/, is introduced which contains only 
soft and hard-collinear fields. While the final effective theory, called SCET//, contains only 
soft and collinear fields. We will then do a two-step matching from QCD — > SCET/ — > 
SCET//. 



A. QCD to SCET/ matching 

In SCET/, the K* meson is taken as a hard-collinear particle and the relevant building 
blocks are Xh c , Xj^, A^ c and h. The velocity of the B meson is defined as v — Pe/m^. The 
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matching from QCD to SCET/ at leading power may be expressed as 

4 4 



Heff - ~^KV» (£ I ds Cf(s)jf(s) +pjdsjdr Cf(s, r)Jf(s, r) 

+ Ids [dr [dtC c {s,r,t)J c {s,r,t)) , (10) 



where C\ A ' B) and C c are Wilson coefficients in the position space. The relevant SCET/ 
operators for the B — > K*£ + £~ decay are constructed by using the building blocks mentioned 
above 2^ |: 

jf = X hc (sn)(l + 75)7^(0) V , J 2 A = X hc (sn)(l + j5 )^L-h(0) , 



a 



J* = X hc (sn)(l + 75)7^(0) , J£ = X hc (sn)(l + 7s)— MO) h»ld , 



J* = X hc (sn)(l+ l5 )4 hc± (rn) h(0) , (11) 

n ■ v 

J 3 B = X hc (sn)(l + n/5h±4hc±(rn)h(0) i'j^i , 

— Tl^ — 

jf = X hc (sn)(l + -f 5 )4hc±(rn) h(0) , 

n ■ v 

J c = X hc (sn)(l + l5 )lx hc (rn) X Ic (an)(l + 7s)^(0) , 

where the operators jf and J? represent the cases that the lepton pair is emitted from 
the b — > s transition currents, while J c represents the diagrams in which the lepton pair is 
emitted from the spectator quark of the B meson. Except the lepton pair, the operators J { ' 



li£ 



it transition currents in SCET, 



23, 



have the same Dirac structures as those of the heavy-to- 

which were first derived in Ref. [la for J A and in Refs. 24 , l25|l for Jf (see, also Refs. 

Ei fir 
1]). In this paper we take the operator basis of |2fiL l3Jj which makes Jf multiplicatively 

renormalized, but we have neglected the operators which contain the Dirac structure J^hc±.l± 

and which do not contribute to the exclusive B meson decays. It is also clear that the 

structure l^y^l arises solely from Q w of the weak effective Hamiltonian. 

Since in practice the matching calculations are done in the momentum space, it is more 

convenient to define the Wilson coefficients in the momentum space by the following Fourier- 
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b 




Ql-B 



(0 

b NQiqX B 

bd M W 

FIG. 1: 0(a s ) contributions to the matching of Qi to A-type SCET currents. The crossed circles 
denote the possible locations from where the virtual photon is emitted and then splits into a lepton 
pair. 

transformations: 

Cf(E,u) = Jds Jdr e l(us+ar)n - p Cf(s,r) , (12) 

C c (E,u) = Jds Jdr J da e l(us+ar)n - p e mn - q C c (s,r,a) , 

with E = n • vn • P/2 and u — l — u. To get the order a s corrections to the decay amplitude, 
we need to calculate the Wilson coefficients Cf to one-loop level and Cf and C c to tree 
level. In the following we will use AjC^ ' to denote the matching results from the weak 

ABC 

effective operators Qj to the SCET currents J { ' ' . With this, the matching coefficients 
from QCD — > SCET/ can be written as 

10 

c <.w> = £ AjC <w, WdiAi)i (13) 
J'=l 

where Hqcd is the matching scale and /x is the renormalization scale in SCET/. 

Each operator of the weak effective Hamiltonian, namely <5i-io> w iU contribute to Cf 
at order a s level, as shown in Fig. 1. But due to the small Wilson coefficients C 3 _ 6 , it is 
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numerically reasonable to neglect the contributions from Qz-§. For the operators Qi,2 and 



Qs, the results can be easily derived from Eqs. (11) and (25) of Ref. 



^1,2^2 (A'qcd) 



qcdJ 



27T 4:71 



(2F 2 (7) + sF? s )C 2 + 2(Ff j + F 2 w /6)d 



,(9) 



(9) 



,(9) 



(2F 2 (7) + F^)C 2 + 2{F^ + if'/fyC! 



^(9) 



qcd; 



A 8 C 2 A (/i 



27T 47T 

(a*qcd) m b (fi QCU ) 
2n An m b 

(a^qcd) m b (fi QCU ) 



(9) 



--(9) 



(9) 



qcd; 



FP + F f 



2FP + f?> 



a 



eff 



(14) 



a 



eff 



2n 47r m fe 

where s = q 2 /ml and m& is the pole mass of the b quark. The current mass m b is related to 
the pole mass at next-to-leading order by 



m b (fi) = m b 



1 + 



a s C F 



(15) 



where C p — 4/3. The functions i^ag! are given in a mixed analytic and numerical form 
in Ref. 32 1. Following the convention of Ref. [14], we also use the "barred" coefficients 
(7j(i=l,...,6) here which are the linear combinations of the Wilson coefficients Cj of the 
weak effective Hamiltonian in Eq. (J3|. The effective Wilson coefficient C| is defined as 
C 8 e// = C 8 + C 3 - C 4 /Q + 20C 5 - 10C 6 . 

For the operators Q7, Qg and Qio, the matchings to the A-type currents give 

a em m b (/i QCD ) 



A 7 Cf = 






2tt 


m b 


a 9 c? = 


2tt 




= 


"em 

2vr 


C3C10 , 



r 1 n e ff 

Og0 7 



AtC-2 4 



A 9 C 2 A 



AioCf 



27r m& 

O^m. / ^ 1 



-ZOio0 7 , 



2vr 

~2tT 



Ca 



c 4 + 



-a a 



eff 



(16) 



-a c 
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To avoid confusion with the Wilson coefficients in Eq. (HI), w e use the notations Cj for the 
matching coefficients, instead of Cj used originally in Ref. (l6j . The explicit expressions of Cj 



up to one-loop order can be read from 
tensor current in 



16 



] . Note that although the operator basis of the 
23j looks slightly different from that of [16l], they are actually the same and 
it is easy to find the relations C$ = C^ 4 and C10 = Cy . The effective Wilson coefficients 
are defined as C 7 e// = C 7 - C 3 /3 - 4C 4 /9 - 20C 5 /3 - 80C 6 /9 and C 9 e// (g 2 ) = C 9 + Y{q 2 ), 
where the function Y (q 2 ) rep resents the contributions of the fermion loops and the explicit 



formula can be found in 




FIG. 2: Tree-level matching of Q, t onto B-type SCET currents. The crossed circles denote the 
possible locations from where the virtual photon is emitted, while the crosses mark the possible 
places where a gluon line may be attached. 



To get the decay amplitude of B — > K*£ + £~ in order a s , the tree-level matching of 
the effective weak Hamiltonian (4) onto B-type SCET currents (11) is already enough, as 
illustrated in Fig. 2. If we use the notation AiqC^ to stand for the matchings of Qi-6 onto 
B-type SCET currents Jf , namely A 16 Cf = Y?j=i A^Cf , we S et from Fi S- 2a that 

A 16 Cf = ( \f^{u, s, ml/ml){C 2 + C,- C 6 ) - V 6 (u, s, 0)C 3 - 

l Fl i(u,s,l)(C 3 + C 4 -C 6 -AC 5 )) , 

a 2 /2 „ 1 „ (17) 

A 16 Cf = ^P— -Fjl («, S, m 2 c /m 2 b )(C 2 + C 4 - C 6 ) - s, 0)(7 3 - 

Z7T TOj, \0 O 

iFjl( M ,S,l)(C 3 + C 4 -a)^) , 
where w is the momentum fraction carried by the strange quark in the K* meson. The 
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functions F 16 '" are defined as 



r , \- i " I I ^-s + 4A V§~ 
r ie (u, s, X) = 1 + — — — r- s j= — arctan ■ 



(l-a)(l-«)\^ ^ VT V-s + AX 

V-l +m - su + AX J\ - (1 - s)u \ , _ . / 2a/T 

arctan — + A L12 



a/1 - (1 - s)w V-l + u- su + AX J \y/T - Vs - 4A 



V§~ + Vs - 4A J sju + a/1 - (1 - s)u - 4A 

-A Li 2 I - 1V^^> || , (1! 

" 1 a/1 - (1 - s)*z - a/1 - (1 - s> - 4A ' 1 



,-,11 , . ,v . 4s / . / V-s + 4A \/.s 
FJg(-u, s, A) = 2s + — — — r (1 — u + ms) = — arctan 



(l-s)(l-u)V I a/F V-s + 4A 



V-l +u - su + AX V 1 - (1 - s)u \ . _ . / 2vT 

arctan = + A Li 2 



l-(l-s)it V-l +M- su + AXJ \ VT - Vs - 4A 



VF + Vs - 4A y " I V 1 - (1 - S ) M + V 1 ~ (1 - - 4A 



A Li 2 V : \ ; = . 19 

" y/l-(l-s)u-y/l-(l-s)u-4\J J 

As a check, it is not difficult to find the following relations 

n ,, A (l-u)S _n , A m 2 .. . , s(l-u)E 

F id u i s i—2) =t±{u,m q ) x — — — , Fh{u,s,—j) =t\\{u,m q ) x — , 

Tn b zivi B m b . . 

where the functions t± t \\ (u, m q ) are defined in Eqs. (27)-(28) in the paper by Beneke et al. [14]. 
We also note that the functions F^(u, s, A) and Ff 6 (u, s, A) are finite as u = 1 — u — > 0, as 
opposed to the functions t± t \\(u,m q ), which are singular as u — > 0. 

Fig. 2d and the operator Q 9 of Fig. 2f, combined with Fig. 2b, will contribute to the 
matching coefficients A^gC^, while the operator Qio of Fig. 2f will contribute to AioC^ 4 : 

A 9 Cf = , (20) 

AioC,f = , A 10 Cf = — ^ — 1-C10 . 

27r mb{l — s) 

Finally, Fig. 2e and Fig. 2c contribute to the matching coefficients 

b = _a g m^ 2(l- M )(l-g) ff B = 

81 2n m 2 b 3s(u + s - us) 8 ' 82 K J 
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We shall now consider the diagrams where the virtual (off-shell) photon is emitted from 
the spectator quark, as shown in Fig. 3. Due to the off-shellness of the quark propagator, it 
is easy to check that Fig. (3d-3f) are of order l/m b suppressed compared with Fig. (3a-3c) 
where the photon is emitted from the spectator quark in the B meson. Therefore at leading 
power in 1/rrib, only the first three diagrams in Fig. 3 are relevant for our analysis. As we 
shall see in the following, all of these three diagrams contribute to the Wilson coefficients of 
the C-type SCET current. 

The annihilation diagram, shown in Fig. 3a, contributes to the matching coefficient C c 
at order for which the calculation is trivial, 

*®C° = \ ("^(Ci + ZC 2 )8 qu + (C 3 + 3<7 4 )) • (22) 

Here q is the flavor of the spectator quark in the B meson and the superscript (0) denotes 
the matching at order a°. At order a s , the diagrams shown in Figs. (3b-3c) also contribute 
to the matching onto the C-type SCET current with the coefficients 

A 8 C C = ° Fas ~ ACl " , 
N c 4n 1 — u + us 

A^C C = 2-f ^ { (C 2 + Ck + C 6 )G(u, s, m 2 Jm 2 b ) + (C 3 + 3C- 4 + 3C 6 )G(u, s, 0) (23) 

+(C 3 + C A + C e )G(u, s, 1) + t(C 3 -C 5 - 15C 6 ) 

12 



where the function G(u, s, A) is defined as 

2 2 m 2 f 1 

G(u, s, A) = - + -In— £ + 4 / dx x(l -x)\n[X-x(l - x)(l -u + us)] . (24) 

3 3 \i l Jo 



B. SCET/ SCET// matching 



As shown in Refs. [2JJ, |25j , which analyzed the form factors in the framework of SCET, one 
may simply define the matrix elements of the A-type SCET/ currents as non-perturbative 
input since the non-factorizable parts of the form factors are all contained in such matrix 
elements. Therefore the explicit matching of to SCET// operators is not necessary here. 

For B-type SCET/ operators, they are matched onto the following SCET// operators 



Of = X c {sn){l + 75 )7^X C (0) Q,(fn)(l - 75 )^(0) 



Of = X c (sn)(l + 75 )-^-|x c (0) Q s (tn)(l + 75 )|?{ s (0) c-^ , 
Of = X c (sn)(l + 7s)7^ X c (0) Q s (tn){l - 75)^(0) h»ld 



(25) 



2 w v /v ""2 
Of = X c (sn)(l + 75 )-^-|x c (0) Q s (tn)(l + 75 )^ s (0) £7^ , 

where we only include the color-singlet operators that have non-zero matrix elements for 
the B — ► K*£ + £~ decay. Again, it is in practice more convenient to do the matching 
calculations in the momentum space, and the Wilson coefficients Df(u,u) can be defined 
by Fourier transforming the corresponding ones Df(s,t) introduced in the position space, 
just like the case in SCET/, 



Df(u,u) = ds dt e- iujn vt e iusn P D?(s,t). (26) 

Following the notations of |26j], the Wilson coefficients Df can be expressed as 

D?(u,u,s,v) = lj^dv Ji ^.ln ^'^ j Cf(v,fi) , (27) 

where the jet functions Ji arise from the SCET/ — ► SCET// matching and it is clear that 
Jx = J 3 = J± and J2 = Ji = J\\- At tree level, using the Fierz transformation in the 
operator basis, 

X C NX S Q S MX C = -ix c (l + 75 )|x c Q S M(1 - i^N% s - ^X c (l - 75 )|x c 

x Q S M(1 + 75 )|iV^ - \x c {l + 7 5 )|7±«X C Q S M(1 + l5 )^NX s , 

(28) 
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D°(u, u, s, /,) = ? ^ ., J C f In !!^LJ) , „ ) C C (E, u, /,) , (32) 



one obtains 

J±( „, „) = Mu , „) = _^_I__,(„ _ „) . ( 29) 
Finally, the C-type SCET/ current is matched onto the SCET77 operator 

O c = X c (sn)(l + 75 )|x c (0) Q 5 (tn)(l + 75)7^1(0) • (30) 

We may similarly define 

D C (iU, U)= Jds Jdt e -^n-vt e iusn-PjjC^^ f y (31) 

with 

D c {u,u,s,fi) = ., , ... 

[yj — q l I — le) \ /j, z 

where e g is the electric charge of the spectator quark in the B meson. At tree level the 

corresponding jet function is trivial, J c = 1. For later convenience, we will define D c = 

D c /(u - q 2 /m b - it). 



C. Matrix elements of SCET operators 

The last step before we can finally get the decay amplitude for the B — > K*£ + £~ decay is 
to take the matrix elements of the relevant SCET operators. For the A-type SCET currents 
(11), one may simply define j3] 

(M(p)\X hc Th\B(v)} = -2E( M (E)tr[M M (n)TM B (v)} , (33) 

where the projection operators are 

M B (v) = -l±i l5 , M K1 (n)=f ± ^, M Kl {n) = -^, (34) 

with being the polarization vector of the K* ± meson. It is then straightforward to get the 
matrix elements of the SCET/ currents Jf as 

_ 77 ^ _ 

{K*t+r\j£\B) = -2E( ± (g? - ie^elJ^e , (K*tr \ 4\B) = -2EQ 



n ■ v 



(K*£ + r\jf\B) = -2EC±(g? - ie^)e\ v h^i , (K*£ + r \ j£\B) = -2E^^—t llxl d , 

(35) 
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where = g^ v — (n fl n u + n fJ- n u )/2 and = e^^VpUa/in • t>). Note that in the above 
equations, we use the convention e 0123 = +1, as adopted in the book by Peskin and Schroeder 



33|. 



For the B-type SCET77 operators (25), although naively the soft and collinear degrees 
of freedom seem to be decoupled, the factorization may be invalidated unless no endpoint 
divergences a ppea r in the convolution integrals Q, Q| The relevant meson LCDAs are 
defined as 



(0\Q a (tn)m a {0)\B(v)) =^M^ / du ; 



tr 



<f>»(u,ii) - -J—itjZiu^) - <ft(u,ii)) TM B (v) 
In - v 



, (36) 



(IT-Cp) |X c (»n)r|x c (0) |0) =Jj^-n-ptr[M K *T} f du e ~<M«,/x) 



where two different ^-distribution amplitudes ((j) K »(u, fi) for T = 1 and <pj^*(u, fi) for 
r = 7j_) with their corresponding decay constants /|-* and f^*(jj), respectively, are involved; 
F(n) is related to the B meson decay constant fg up to higher orders in by (3^ 

/bv* = F W (l + ^> (3 In ^ - 2)) . (37) 

With the above LCDAs, the matrix elements of the operators Of can be written as 

(K*m- \C?0* \B) = - F MfJ? {1 _ §){gT _ ie^elJ^i !°° 



x / du f K * {lA<t>K* (u, n) [ dvJ±(u,v,\n mbUJ (\ — — , /i)Cf (v, /x) 
Jo V> 

4 



-(1 - S)G/T - 0+ ® /ki^! ® Jl ® Cf 



(K*£ + r\CiO B \B) = ^^-(1 - s) h,£ 0* ® f Kl <f> Kl ®J ]{ ® Cf , 

(38) 

while for the matrix element of CfO^(CfOf), it can be obtained by simply replacing the 
lepton current ij^i on the right hand side of the above equations by £7^75^ and also replacing 
c? - Cf {C* -> Cf). 

The matrix element of 0° is obtained likewise, with the result 

(K*en~ \D c O c \B) = - F ^\ m B 2 {1 _ s)^h/ ^ ® /«&q ® 5 C . (39) 
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Since (j>^(u) does not vanish as uj approaches zero, the integral f du (j) B {uj)/{u) — q 2 /rrib) would 
be divergent if q 2 —>■ 0. This endpoint singularity will violate the SCET// factorization, that 
is why we should restrict our attention to the kinematic region where the invariant mass of 



D. Resummation of logarithms in SCET 

In the above analysis a two-step matching procedure QCD — > SCET/ — > SCET// has been 
implemented. This introduces two matching scales, \Xh ~ at which QCD is matched onto 
SCET/ and fii ~ y/m^K at which SCET/ is matched onto SCET//. Thus, with the SCET/ 
matching coefficients at scale fih, one may use the renormalization-group equations (RGE) of 
SCET/ to evolve them down to scale \xi and then match onto SCET//. The large logarithms 
due to different scales are resummed during this procedure. Note that the meson LCDAs 
may be given at another scale fii, and, in principle, one should also use the RGE of SCET// 
to run the corresponding matching coefficients from fii down to fii- But since in B decays 
the scale /i; ~ 1.5GeV is already quite low, we may just take the meson LCDAs at the 
scale Hi in this paper for simplicity and thereby avoid the running of the SCET// matching 
coefficients. 

Furthermore, one should note that for the A-type SCET currents, only the scale fih is 
involved since it is not necessary to do the second step matching of SCET/ — > SCET//. 
Similarly, we may choose the nonperturbative form factors (jii at the scale \ih and avoid 
the RGE running of the A-type SCET/ matching coefficients. For the B-type currents, the 
RGE of SCET/ can be obtained by calculating the anomalous dimensions of the relevant 
SCET operators, which has been done in 3lj . where the matching coefficients at any scale 
fi can be obtained by an evolution from the matching scale fih as follows 



with the subscript T =_L, || and the functions a(/i/ l ,/i) and S(fj,h,fi) are given in Eq. (66) of 
Ref. 0|. Note that in the above equation one should use the subscript T =_L for j = 1,3, 



the lepton pair is not too small, say q 



■ 2 > 1 GeV 2 . 




(40) 
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while T =|| for j = 2,4. The evolution kernel U[(E , u, /x^, //) obeys 

^rfh/T'^ = / dl/ ^ r(2/ ' U)d3r{E ' V ' ^ ^ + U( - U W(E, u, A* fc , /*) , (41) 
with the initial condition U[(E, u, fi^, /i/J = Cf(E, u, /z^). Again, the functions Vr(y, u) and 



u;(it) are defined in |31j. In the next section on phenomenological application, we will solve 
the above integro-differential equation numerically. 

Finally, for the C-type SCET current J c , its anomalous dimension just equals the sum 
of the anomalous dimensions of the K* meson LCDA (fix* and the B meson LCDA <p B . 
However, as the evolution equation of 0^ is still unknown, we will not resum the pertur- 
bative logarithms for the J c current in this paper. Numerically the contribution from the 
J c current to the decay amplitude is small. Furthermore, as we will see later, the J° cur- 
rent is completely irrelevant for the forward-backward asymmetry of the charged leptons. 
Therefore, this treatment has only minor impact on our phenomenological discussion. 

III. NUMERICAL ANALYSIS OF B -» K*£+£- 

We are now in the position to write the decay amplitude of B — > K*£ + £~, using the 
similar notations adopted in |l^ |. 

d 2 T 



G 2 F \v t :v tb \ 2 faemy,, , 

dq 2 d COS 9 " 1287T 3 WW b a k* \ m y 



2 x 



2 



2Ci(l + cos^)^(|^| 2 + (C i L ) 2 ) (42) 

m B 

_8 C icos^i?e(^)C 1 J + Cjf(l " cos 2 ^)(|C|| 2 + (Cf ) 2 ; 
m B 

with meA^*/2 being the 3-momentum of the K* meson in the rest frame of the B meson, 



i-iV-2^ fi+4w^ 



2 I 9 I 2 1' 4 

171 B J m B V m B J m B 



1 1/2 

(43) 



The angle 9 denotes the angle between the momenta of the positively charged lepton and the 
B meson in the rest frame of the lepton pair. Note that in the above equations the leptons 
are taken in the massless limit and the K* meson mass is kept nonzero only for A^* , which 
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arises from the phase space. The "effective" Wilson coefficients C Q and C 10 are given by 
C 1 - = — (c A + mB fs ^ ® ® J± ® ^ 



2tt / , m B / B $! ® 4.^. (8) J,, <g> 
c 9 — - o 2 + 



q 2 f B u<P^/(u-q 2 /m b -it)® f\,<p\,®D c ' 



^era 



4m B C|| 

.1 



(44) 



n \\ 2tt / A m B futf ® fk.fik* ®J\\ ®Cf 
(-10 — - W + 



where Cj A ' B and Z) c are defined in Eqs. (13) and (32), respectively. The above expressions 
are valid at leading power in 1/mb and to all orders in a s . But in this paper we only calculate 
explicitly the "effective Wilson coefficients" at one- loop order. At this order our results are 
quite similar to those of ^| using the large-energy limit of QCD. The main phenomenological 
improvement is that for the hard scattering part, the matching coefficients Cf are evolved 



from the scale /^ ~ C?(m&) down to ~ y/m b A h , during which the perturbative logarithms 
are summed. Here, represents a typical hadronic scale. Note also that the definitions of 
the soft form factors £m| in SCET are different from those of Ref Q|, therefore the explicit 
expressions for Cf are also different from the coefficients C®' 1 appearing in which are 
related to the form factor corrections. 

In terms of the helicity amplitudes for the decay B — > K*(— > K + 7r)£ + £~, the double 
differential distribution d 2 B/d cos 8 + ds is given in Eq.(44) of Ref. [36]. This requires the 
helicity amplitudes, \H {s)\ 2 = \H^{s)\ 2 + \H^(s)\ 2 , \Hb R (s)\ 2 and \H^ R (s)\ 2 . While the 
amplitudes H+ R (s) are both power suppressed in l/m b and numerically small, the expres- 
sions for the others in SCET are given below: 

w = ^(i-^-) 2 (rii 2 +(cj ) 2 Kf , 



B 



(45) 

\H-*\ 2 =t£-\?\#±c&\ 2 <i- 

m B 

Note that the dependence on the soft form factors factorizes in Q and (]_ for the helicity 
components \H \ 2 and \H^ ,R \ 2 , respectively. Since a similar analysis in terms of the helicity 
amplitudes of the charged current decay B — ► p(—> itit)£ + i>i can be performed, the ratios 
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Ro(s) and R-{s) of the two differential distributions (in B — > K*(— > Kn)£ + £~ and S 
p(— >• 7T7r)£ + z/£) have lot less hadronic uncertainties, as these ratios (see Eq. (76) in Ref. 
for their definition) involve estimates of the SU(3)-breaking in the soft form factors. The 
point is that the ratios / and (*[ / are more reliably calculable than the form factors 
themselves. 



A. Input parameters 



To get the differential distributions numerically, some input parameters have to be spec- 
ified. For the calculation of the Wilson coefficients, the relevant parameters are chosen as 



M w = 80.425 GeV , sin 2 9 W = 0.2312 , A. 



(5)_ 
MS 



217l| MeV 



and m? t ole = (172.7 ± 2.9) GeV, updated recently by the Tevatron electroweak group 



(46) 

Q 



Numerical values of the Wilson coefficients, evaluated at scale \i = rrib = 4.8 GeV, with the 
three-loop running of a s and the input parameters fixed at their central values given above 
are shown in Table I. Note that the NNLL formula for Cg can be found, for example, in 



the appendix of |l4j, while the re 
have been calculated recently in 



evant elements of three-loop anomalous dimension matrix 



39 



40] 



TABLE I: The leading-logarithmic (LL) and next-to-leading-logarithmic (NLL) Wilson coefficients 
evaluated at the scale mj = 4.8 GeV. For Cg^o, they are also given in the NNLL order. 





LL 


NLL 




LL 


NLL 


NNLL 


Gi 


-0.2501 


-0.1459 


G 6 


-0.0316 


-0.0388 






1.1082 


1.0561 


C e 7 " 


-0.3145 


-0.3054 




c 3 


0.0112 


0.0116 


r eff 


-0.1491 


-0.1678 




G 4 


-0.0257 


-0.0337 


Cg 


1.9919 


4.1777 


4.2120 


G 5 


0.0075 


0.0097 


Gio 





-4.5415 


-4.1958 



The CKM factor \Vu\Q,\ ~ (1 - \ 2 /2)\V cb \ is estimated to be 0.0403 ± 0.0020 by taking 
\Vcb\ = 0.0413±0.0021 |4l| and A = 0.2226. For the B meson lifetimes, we use t b + = 1.643 ps 



and tbo = 1.528 ps |4l|. The pole mass mj, is chosen to be 4.8 GeV. The ratio of the charm 
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quark mass over the b-quark mass is taken to be m c /rrib = 0.29 ± 0.02. For the matching 
scale from SCETV to SCETYr, we use fii = \fm^K} l ~ 1.5 GeV. 

The hadronic parameters for the decay B — > K*£ + £~ include decay constants, light-cone 
distribution amplitudes (LCDAs) and the soft form factors. The B meson decay constant can 
be estimated by QCD sum rules or lattice calculations, here we take Jb = (200 ± 30) MeV. 
For the K* meson, experimental measurements give 37] f\« = (217 ± 5) MeV while the 
most recent light-cone sum rules (LCSRs) estimate 42^ is /^*(1 GeV) = (185 ± 10) MeV. 
Note that f^, obeys the scale evolution equation /^»(/i) = /^.(/io)( a s(/ i )/ a s(/ i o)) 4/ ' 23 - 

The B meson LCDAs enter into the decay amplitudes only in terms of the integrated 
quantities and X B 1 _(q 2 ) defined as by the following integrals 



'dou 



Mi 



10) 



0f(w) 



du- 

lu - q 2 /m, b - ie 



(47) 



Therefore, it is not necessary to know the details about the shape of <f)^_(uj). The most recent 
estimate gives Q = (1.86 ±0.34) GeV" 1 at the scale \i = 1.5 GeV. However, \ B ]_{(f ) 
does require the knowledge of <p^_(uj), about which we know very little. Fortunately, X B 1 _(q 2 ) 
only appears in the annihilation term which plays numerically a minor role in the B — > 
K*£ + £~ decay. To be definite, we adopt a simple model function 3^| ^(co) 
with Un 1 ~ 3 GeV" 1 . 



The K* meson LCDAs may be expanded in terms of Gegenbauer polynomials: 



0^'i'(u,/i) = 6ii(l — u) 



1 



n=l 



(48) 



However, the coefficients a n are largely unknown. Following [44], we shall ignore the terms 
a^'" (n > 2). For ai^, we omit their scale dependence and estimate in a conservative manner: 



a 



0.1 ± 0.1, a 



0.1 ± 0.1. We note that recently the first Gegenbauer moment of 



the K* meson has been revisited in LCSRs 



42| which gives smaller uncertainties. 



There are only two independent B — > K* form factors in SCET, namely (±(q 2 ) and (\\(q 2 )- 
They are related to the full QCD form factors as discussed in |3l|. The current knowledge 
of these form factors is fragmentary. For instance, (± may be extracted from V B ^ K * 



a(g 2 ) 



q(o) 



1 — q 2 /rriy 



+ 



1 



q 2 /m 2 



(49) 



"Vfit i 

with r\ = 0.923, r\ = -0.511, m v = 5.32 GeV and m 2 Vfit = 49.40 GeV 2 . Note that the 
g 2 -dependence above is the same as that of V B ^ K * (q 2 ), calculated in H using LCSRs. 
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However, analyses of the radiative B decays and 
the semi-leptonic B decay B — > plv 47] imply that the LCSRs overestimate the B — > V 
form factors si gnifi cantly. We use the radiative B — > K*^ decay, which has been measured 
quite precisely |4l|]: B(B° — > K*°^y) = (4.01 ±0.20) x 10~ 5 , to normalize the soft form factor 



at q 2 = 0. In SCET, it is straightforward to get the decay amplitude of B — > K*^ from the 
B — > K\i + i~ decay, by taking the limit q 2 — > 0. Then, using the input parameters from 
Table II, we obtain C±(0) = 0.32 ±0.02. Here the error is mainly from the CKM factor V ts V t l 
and the experimental uncertainty of the branching ratio B(B° — > K* 0r y). This estimate is 



consistent with the result of Ref. 



261 ] . but significantly smaller than the number 0.40 ±0.04 



we get from LCSRs. In our numerical analysis, we will choose the value (±(0) = 0.32 ± 0.02 
determined from the radiative B decays, but assume that the g 2 -dependence of C-i_(<? 2 ) can 
be reliably obtained from the LCSRs. 

For the longitudinal soft form factor £m , unfortunately there is no quantitative determina- 
tion from the existing experiments, though this may change in the future with good quality 
data available on the decay B — > piuf. Using helicity analysis, one can extract (^(q 2 )] com- 
bined with estimates of the SU(3)-breaking one may determine C|f"*(<? 2 )- Not having this 
experimental information at hand, one may extract C||(<? 2 ) from the full QCD form factor 

a s (m b )C F 



ab^k 



-4, 



B^K* 



(a 2 



47T 



21n 2 [l 



9 

-ln[l 

s 



s] ±2 Li 2 [s] 



4±^ 
12 



4(1 



BS 



i,B 



II Jl 



J\\ ® [ — ] e s( ^ J dy U\\ (v, y, p h , W ) 



(50) 



with s = q 2 /m 2 B . LCSRs estimate |4J A^ K *(0) = 0.374 ± 0.043 with the g 2 -dependence 

1.364 0.990 



a b->k* 



(q 2 



(51) 



1 - q 2 /m 2 B 1 - g 2 /36.78GeV 2 ' 
From which we get C||(0) = 0.40 ± 0.05, using the input parameters discussed above and/or 
listed in Table II. Its g 2 -dependence is drawn in Fig. 4. 

Alternatively, Gil 2 ) may also be determined from the following relation, 



Em B {V - A 2 



>B^K* 



m K *(m B + m K *) 

1 - 2s 



4(1 -s) 



f. 



BS 



" (21n 2 [l - s] - 2 Ml -s}+2 Li 2 [s] + 6 + ^)] C||(? 2 

— e ^'«) / dyU^v^fih,^) . 

/ Jo 

(52) 
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TABLE II: Numerical values of the input parameters and their uncertainties used in the phe- 
nomenological study. 



M w 


80.425 GeV 


sin 2 0w 


0.2312 


pole 
m t 


(172.7 ± 2.9) GeV 


A^L 

MS 


(217±2J) MeV 


\ V tsVl\ 


(40.3 ± 2.0) x 10~ 3 


aemimb) 


1/133 


TUB 


5.279 GeV 


pole 
m h 

O 


4.8 GeV 


T B+ 


1.643 ps 


T B 


1.528 ps 


m c /m b 


0.29 ±0.02 




1.5 GeV 


X B ] + (1.5 GeV) 


(1.86 ±0.34) GeV" 1 


fB 


(200 ± 30) MeV 


Cl(0) 


0.32 ±0.02 


Cn(o) 


0.40 ± 0.05 


fk*0- GeV ) 


(185 ± 10) MeV 


/i 


(217 ± 5) MeV 




0.1 ±0.1 


a 2 


0.1 ±0.1 



With the input V B ^ K * (0)-A$^ K * (0) = 0.152±0.057 from LCSRs, we obtain Cy (0) = 0.42± 
0.16, which agrees with the range extracted from A B ^ K \ We will use C||(0) = 0.40 ± 0.05, 
obtained from its relation to the full form factor Aq~* and the LCSR, as discussed above. 
Fig. 4 shows the g 2 -dependence of both soft form factors C±,||(? 2 )- However, since the analysis 
of the semileptonic decay B — > plv 4?J suggests that both the transverse and longitudinal 
form factors might be overestimated by LCSRs, we will also consider, as an illustration of the 
non-perturbative uncertainties, the value C||(0) = Cl(O) = 0.32 with all the other parameters 
taken at their central values. 



B. Numerical solution of the SCET/ evolution functions 

As we discussed in Sect. II. D, the B-type matching coefficients Cf should be run from 
the scale = 4.8 GeV down to fii = 1.5 GeV, with the evolution kernel Ur(E, u, fih, n) 
obeying the integro-differential equation (}4l"j) . To solve this equation numerically, it is more 
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12345678 

FIG. 4: The (^-dependence of the soft form factors C±,\\(q 2 )- The solid curve represents Ci(<7 2 ), 
while the dashed curve represents C\\(q 2 )- We have rescaled the transverse form factor at q 2 = 0, 
to be consistent with the experimental measurements of the B — ► K*j decay rate. 




g 2 [GeV 2 ] 



convenient to define the following evolution functions, 



Uy(E,u,fi h ,n)= / dvUr(u,v,n h ,n) , 
Jo 

u + s — us f 1 1 — V , . 

Uy{E,u,fj lh ,fj,) = — / dvUr{u,v,/j, h ,/j,) — — -, (53) 

1 — u J v + s — vs 

l4 c) {E,u,fJL h ,fj,) = / dvU r (u,v,n h ,n)- 1(1 



F[ 6 (u,s,m 2 Jm 2 ) ' 

where T =_L, || and the functions F{g '(it, s, m 2 /m 2 ) are denned in Eqs. (fTHj) and (JTHJ). Note 
that at the quark level, the K* meson energy is related to s by E = m^l — s)/2 in the rest 
frame of the 6-quark. With such definitions, the above evolution functions are normalized 
to one at the scale u^: U$?' b ' c ' (E,u, [ih, Uh) — 1> an d the QCD parameter A§= would be 
the only input for their numerical evaluations. The matching coefficients Cf at scale \i% can 
then be written as 

— e^^U^^E^^^^iCfiE^iih) , (54) 

where we should use the superscript (a) for A 7 9iW C^, the superscript (6) for A%Cf and the 
superscript (c) for A^Cf. For the subscript T, one should use T =_L for j — 1, 3 and T =|| 
for j = 2,4, which is the same as the convention of Eq. (j40|) . Note that for the evolution 
of A^Cf , we have taken into account the fact that the term F[ 6 (u, s, m 2 /m\) is dominant 
due to the large Wilson coefficient 
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It is then straightforward to get the following evolution equations 
dU^ a) (E,u,fi h ,fi) ''' 



dln/i 



/ dy yV r (y,u)U^ a \E,y,fi h ,fi) + u(u)U^ a) (E,u,fi h ,fi) 
Jo 



— = / dyyV r {y,u)- ( —Uy(E,y,fj, h ,fi) + 

din// J (1 -u)(y+ (1 -y)s) 

uj(u)U^ b) (E,u, ji h , ji) , ^ 



dU^ c) {E,u,n h ,n) f 1 F[ 6 (y,s,m 2 c /m 2 b )~ (c) , 

Lo(u)U^ c) (E,u,fi h ,fi) . 

To get the numerical solutions of the above integro-differential equations, we will perform 
the scale evolution in one hundred discrete steps. While from the scale /i n to /t n +i, the 
convolution integral is evaluated for three hundred different values and discrete s values of 
5s = 0.01 in the interval s £ [0.04,0.35]. The function Ur(E, u, /t^, // n +i) is obtained from 



(5) - 

a fit to these values. Taking = 217 MeV, the numerical results of these evolution 

r( a ), 



functions are shown in Fig. 5. Note that the function ' (E,u, fj,h, /*) actually does not 
depend on the energy E, as shown in Fig. (5a). In fact, it is just the same function as 
Ur(u, jJhi defined in Eq. (5.23) by Neubert et al. [3Jj. The function C/ij is not shown 
in Fig. 5, since it does not enter into the decay amplitude at the one-loop level, due to 
A 8 Qf = 0. While for the complex functions only the absolute values of the functions 
are plotted. 

C. The dilepton invariant mass spectrum and the forward- backward asymmetry 

Experimentally, the dilepton invariant mass spectrum and the forward-backward (FB) 
asymmetry are the observables of principal interest. Their theoretical expressions in SCET 
can be easily derived from Eq. 



dBr _ G 2 F \V*V tb \ 2 , aem y 3 q 2 , 

{yCl^d^l 2 + (O 2 ) + 4 -( 2 (\cl\ 2 + (Cf ) 2 )} 



(56) 
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(a) 



(b) 




(c) (d) 

FIG. 5: Numerical values of the functions Uj?' b ' c \E,u, Hh, Hi), evolved from Hh = 4.8 GeV down 
to Hi = 1.5 GeV, the relevant parameters are taken at their central values. For the upper-left plot, 
the solid line denotes f/j°' while the dashed line denotes . For the lower plots, since f/p are 
complex functions, we only show their absolute values. 



dA FB = i / r 1 ^^ d 2 v _ r° dcQs9 tPr \ 

dq 2 dT/dq 2 \J dq 2 dcos9 J _ x dq 2 d cos 9 J 

= -6(q 2 /m 2 B )( 2 ± Re(Ci)Ct (57) 

<?/rn%)&{\<£\ 2 + {&) +C|f(l4'l 2 + (C} ) 2 ) ' 
With the input parameters listed in Table II, the decay spectrum and the FB asymmetry 

are shown in Fig. 6 and Fig. 7, respectively. In our calculation we have dropped the small 

isospin-breaking effects, which come from the annihilation diagrams, and take the spectator 

quark as the down quark in Eqs. To estimate the residual scale dependence, we 

vary the QCD matching scale Hh by a factor v^2 around the default value [ih — m b- Note 

that the soft form factors Ci,]|((Z 2 ) defined in SCET are actually scale dependent, which effect 
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FIG. 6: The differential branching ratio <i8(5° -> K*°e + l~)/dq 2 in the range 1 GeV 2 < <? 2 < 
8 GeV 2 . In the left plot, the solid line denotes the theoretical prediction with the input parameters 
taken at their central values, while the gray area between two dashed lines reflects the uncertainties 
from input parameters and scale dependence. In the right plot, the soft form factors are normalized 
as C||(0) = Cl(0) = 0.32, while all the other parameters are chosen at their central values. 
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FIG. 7: The differential spectrum of the forward-backward asymmetry dAps(B — > K*£ + £~)/dq 2 
in the range 1 GeV 2 < q 2 < 8 GeV 2 . Here the solid line denotes the theoretical prediction with 
the input parameters taken at their central values, while the gray band between two dashed lines 
reflects the uncertainties from input parameters and scale dependence. The dotted line represents 
the LO predictions, obtained by dropping the 0(a s ) corrections. 



has been taken into account in our error analysis. 

Restricting to the integrated branching ratio of B — > K*t 



in the range 1 GeV < q 2 < 
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7 GeV 2 , where the SCET method should work, we obtain 

7 GeV 2 

2 dBr{B — » K* £ £ ) _ 0O+ q,57i +0.301 +o.is\ v in -7 /con 

a Q ^ - l 2 - y2 -0.50K|| -0.28ICKM _ .20J X iU ■ V°°) 

1 GeV 2 

Here we have isolated the uncertainties from the soft form factor £y and the CKM factor 
|VkSVt&|. The last error reflects the uncertainty due to the variation of the other input 
parameters and the residual scale dependence. If the smaller value for the longitudinal form 
factor C||(0) = 0.32 is used, as shown in Fig. (6b), the central value of the branching ratio 
is reduced to 2.11 x 1CT 7 . For B° decay, the branching ratio is about 7% lower due to the 
lifetime difference: 

7 GeV 2 

dBr(B° — > K*°£ + £~) 
2 aniyn — » i\ 11 ) +(X53| +0 28 | +0 .17\ v 1fr -7 / kq\ 

aq ^ ~~ l 2 -' 2 -0.47K|| -0.26ICKM _ .19J X iU ■ l^J 

1 GeV 2 

To compare with the current experimental observations, it was proposed in Ref. ^| to 
consider the integrated branching ratio over the range 4 GeV 2 < q 2 < 6 GeV 2 , for which we 
get (0.92+°; 2 g) x 1(T 7 . This is smaller than the number (1.2 ± 0.4) x 10" 7 obtained in Ref. 
|l4j, which is mainly due to the fact that the most recent LCSRs estimation |4jJ prefers the 
form factor A^~* K * to be smaller. Experimentally one of the Belle observations of our 
interest is 

8 GeV 2 

/dBr(B — > K*£ + £~) 
dq 2 " = ( 4 - 8 -L2| S tat. ± 0.3| syst . ± 0.3| mode i) x 10~ 7 , (60) 

4 GeV 2 

for which we predict (1.94lg 40 ) x 10 -7 . This is smaller than the published Belle data by 
a factor of about 2.5. But at this stage, it is still too early to conclude that one should 
change some theoretical input significantly to be consistent with the experimental data. 
For instance, the BaBar collaboration measures the total branching ratio of B — > K*£ + £~ 
o be P| (7.8]ti'7 ± 1.2) x 10~ 7 , which is about twice smaller than the Belle observation 
"\ (16. 5+2:2 ± 0.9 ± 0.4) x 10 -7 . This implies that, if finally the total branching ratio of 
B — ► K*£ + £~ is found to be closer to the BaBar result, the partially integrated branching 
ratio over the range 4 GeV 2 < q 2 < 8 GeV 2 could be lowered to around 2.3 x 10~ 7 , which 



is consistent with our estimate (1.941q; 40 ) x 10 7 within the stated errors. We look forward 
to experimental analyses from BaBar and Belle based on their high statistic data. 
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One of the most interesting observables in the decay B — > K*£ + £~ is the location, 
where the FB asymmetry vanishes. It was first noticed in the context of form factor models 
in j3] and later demonstrated in Q], using the symmetries of the effective theory in the 
large-energy limit, that the value of q$ is almost free of hadronic uncertainties at leading 
order. From Eq. (J57)) . it is easy to see that the location of the vanishing FB asymmetry is 
determined by Re(Cg) =0. At the leading order, this leads to the equation C$ + + 
Re(Y(qD) = 0. Including the order a s corrections, our analysis estimates the zero-point of 
the FB asymmetry to be 

ql = (4.07i» GeV 2 , (61) 
of which the scale-related uncertainty is A(gQ) sca i e =1q;o5 GeV 2 for the range m&/2 < fi h < 



2rrib together with the jet function scale /i/ = yj \ih x 0.5 GeV, as used in the paper by Beneke 
et al. |l4|. Since no reliable estimates of the power corrections in 1/mj, are available, we 
should compare our results with the one given in Eq. (74) of Q], also obtained in the absence 
of l/rrib corrections: q$ = (4.39^35) GeV 2 . Of this the largest single uncertainty (about 
±0.25 GeV 2 ) is attributed to the scale dependence. While our central value for q$ is similar 
to theirs, with the differences reflecting the different input values, the scale dependence in 
our analysis is significantly smaller than that of ^| . This improved theoretical precision on 
q\ requires a detailed discussion to which we now concentrate in the rest of this section. 

As already stated in the introduction, the expressions for the differential distributions in 
the decay B — > K*£ + £~ derived here and in Q| are similar except for the definitions of the 
soft form factors and the additional step of the SCET logarithmic resummation incorporated 
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in our paper. This resummation has also been derived in the existing literature 
However, its effect on the scale-dependence of q$ has not been studied in sufficient detail. 
With the SCET form factors (±(q 2 , fi) and C||(<? 2 > A*) defined in Eq. (33) here, which are scale- 
dependent quantities, and neglecting the resummation effects consistently in both the decays 
B — > K*£ + £~ and B — > K*j, the scale uncertainty is increased, with q% = 4. 12+^07 GeV 2 . 
We draw two inferences from this numerical study: (i) Incorporating the SCET logarithmic 
resummation helps in the reduction of scale dependence in (?q, (ii) A(gQ) scale = lo!o7> obtained 
by dropping the resummation effects is still significantly smaller (by a factor 2) compared 
to the corresponding uncertainty A(gQ) sca i e = ± 0.25 GeV 2 calculated in Ref. This 
deifference, as argued below, is to be traced back to the different definitions of the soft 
form factors used by us for the SCET currents and the corresponding quantities employed 
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by Beneke et al. Q| in the QCD factorization approach. The results in Ref. j3| are, 
however, formally equivalent to the so-called "physical form factor" (PFF) scheme in SCET, 
as discussed subsequently by Beneke and Yang [49]. Thus, the scale dependence of the 
distributions in B — > K*£ + £~, in particular of g^, is related also to the definitions (or 
scheme dependence) of the form factors in effective theories. The PFF-scheme is one such 
choice, but this choice is by no means unique. 

Concentrating on the transverse form factor, relevant for of the FB asymmetry, in the 
PFF scheme, the corresponding SCET/ form factor (where we have now added a suffix 
P for this scheme) is defined as 

CI = m * V , (62) 
m B + m K * 

where V is one of the physical form factors in the decay B — > K*£ + £~ in full QCD. In 
contrast, in our paper, the soft SCET form factors are defined in Eq. (JHEJ). These two defi- 
nitions can be related to each other by = C 3 (j_, where the expression for the perturbative 
QCD coefficient C 3 is given below (C 3 is called Cy 40 " 11 in Since the decay amplitude 

should be independent on how one defines the soft form factors, one must have 

Ci P tl = Ci(± =► Ci p = Ci/C, , (63) 

Since C 3 = 1 + 0(a s ), by expanding Cg /C 3 to order a s , one obtains 
C 1 - 



'9 



i-d-c 3 



2VT ( c a | QW (1 C 3 )(-C e 7 ff I C° ff ) I m B fB<t>l®fK^®J±®Cf 



41n - 4 + ^ ^ ln(l - s) 

jJL 2 S 



(64) 



a em V 1 2tt v ' ' I Cj 

= ci ff + ?C 7 e// (l + %^ 

which agrees with the expression for C^" P in Eq.(40) of Q| (called Cg t ±(q 2 ) there). We recall 
that to determine qfi, we solve the equation Re Cg = 0, where now the quantity Cg is defined 
as follows 

Ci = CMC e g ff + lc e 7 ff ^CM + ... , (65) 
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with the QCD coefficients H (Cg is called C ( T A0)2 in Q 



+2 ln 2 (l - s) + 2Li 2 (s) + Y2 + (| - 3 ) ln (l - s) + 6 



5 9 = 

4?r 



2\ n z ( A] - (41n(l - i) - 7) In ( — 
m 6 y V m & 



7T 2 

+2 ln 2 (l -a) -2 ln(l - s) + 2Li 2 (s) + — + 6 



(66) 



The ellipses above denote the terms which are the same for Cg P and Cg. The functions 
multiplying the effective Wilson coefficients Cg and appearing in Cg P and Cg in 

Eqs. (64) and (65), respectively, lead to different scale-dependence for q^. 

Our result for using the SCET form factors has been given above in Eq. (61) with 
the scale-dependent uncertainty A(g 2 ) sca i c =^05 GeV 2 . Note that we have considered in 
a correlated way the scale-dependence of (±([i,q 2 ) in our analysis. To illustrate this, we 
use the experimental data on the branching ratio of B — > K*j and the central values of 
the other input parameters given in Table II, which yields the following scale-dependence 
of the relevant form factor: £j_(0,/i = 1m^) = 0.34 and £j_(0, /z = = 0.30. In solving 

the equation .RefCg 1 ] = 0, relevant for the zero-point of the FB asymmetry in the decay 
B — ► K*£ + £~, we have factored in the scale-dependence of (±(fi,q 2 )- We do a similar 
numerical analysis of q^ in the PFF-scheme, where the corresponding form factor (±(q 2 ) 
is scale-independent, and incorporate the effect of the logarithmic resummation in both 
the B — > K*7 and B — > K*£ + £~~ decays. Solving now the equation Re[Cg P ] = 0, using 
the central value of the soft form factor ( p (0) obtained from the analysis of the B —>■ K*^ 
branching ratio: ( p (0) = 0.28, and with all the other parameters fixed at their central values 
given in Table II, we find that in the PFF-scheme g 2 = 3.98 ± 0.18 GeV 2 . Had we dropped 
the resummation effect, we would get q% = 4.03 ± 0.22 GeV 2 , where the scale uncertainty 
A(gg) sca ie = ±0.22 GeV 2 , derived here in the PFF-scheme, is consistent with the number 
A((?o)scaie = ±0.25 GeV 2 obtained in jl^. Therefore, we conclude that the difference in the 
estimates of the scale dependence of g 2 here and in Ref. is both due to the incorporation 
of the SCET logarithmic resummation and the different (scheme-dependent) definitions of 
the effective form factors for the SCET currents and the ones used by Beneke et al. [l^ . 
Using the SCET form factors defined in Eq. (J33j) in this paper, we find that the scale-related 
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uncertainty A(gQ) sca i e is reduced than in the PFF-scheme of Beneke et al. jl4j. One expects 
that such scheme-dependent differences will become less marked after incorporating the 
0(a 2 ) effects in the decay distributions for B — > K*£ + £~ . Our comparative analysis hints 
at rather large 0(a 2 ) corrections to q 2 in the PFF-scheme and a moderate correction in the 
SCET analysis carried out by us in this paper. Since the value of q 2 offers a precision test of 
the SM, and by that token provides a window on the possible beyond-the-SM physics effects, 
it is mandatory to undertake an 0(a 2 ) improvement of the current theory of B — * K*£ + £~ 
decay. As power corrections in \jm b have not been considered here, although they are 
probably comparable to the 0(a s ) corrections as argued in a model-dependent estimate 
of the 1/rrif, corrections by Beneke et al. Q|, it also remains to be seen how a model- 
independent calculation of the same effect the numerical value of q 2 . 

IV. SUMMARY 

In this paper, we have examined the rare B decay channel B — > K*£ + £~ in the framework 
of SCET, where the factorization formula holds to all orders in a s and leading order in 
Making use of the existing literature, we work with the relevant effective operators in SCET 
and the corresponding matching procedures are discussed in detail. The logarithms related 
to the different scales Hh = rn b and fii = \/m b A h are resummed by solving numerically 
the renormalization group equation in SCET. We then give explicit expressions for the 
differential distributions in q 2 for the decay B — > K*£ + £~ including the 0(a s ) corrections. 
In the phenomenological analysis, we first discuss the input parameters, especially how to 
extract the soft form factors C±,||(? 2 ) from the full QCD form factors and also the constraints 
on C_i_(0) from the experimental data on the B — > K*^ decay. Using the g 2 -dependence of the 
form factors from the LCSRs and the normalization (±(0) = 0.32 ± 0.02 and C||(0) = 0.40 ± 
0.05, we work out the differential branching ratio and the forward-backward asymmetry as 
a function of the dilepton invariant mass. In the region 1 GeV 2 < q 2 < 7 GeV 2 , where the 
perturbative method should be reliable, our analysis yields 



1 GeV 2 

which can be compared with the B factory measurements in the near future. The largest 
uncertainty in the branching ratio is due to the imprecise knowledge of C||(9 2 )- We have 



7 GeV 2 




(2.92l°f ) x 10 



,-7 



(67) 
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illustrated this by using a value (\\ (0) = 0.32, which reduces the central value of the branching 
ratio to 2.11 x 10~ 7 . We point out that precisely measured ^-distributions in B — > K*£ + £~ 
and B — > ptvn would greatly reduce the form-factor related uncertainties in the differential 
branching ratios. The FBA is less dependent on the soft form factors, and the residual 
parametric dependencies are worked out. We estimate the zero-point of the FBA to be 
(/q = (4.07^0^3) GeV 2 . The stability of this result against 0(a 2 s ) and 1/mj, corrections 
should be investigated in the future. 
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